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Abstract. The origin of brown dwarfs (BDs) is an important component of the 
theory of star formation, because BDs are approximately as numerous as solar 
mass stars. It has been suggested that BDs originate from the gravitational frag- 
mentation of protostellar disks, a very different mechanism from the formation of 
hydrogen burning stars. We propose that BDs are instead formed by the process 
of turbulent fragmentation, like more massive stars. In numerical simulations of 
turbulence and star formation we find that gravitationally unstable density peaks 
of BD mass are commonly formed by the turbulent flow. These density peaks 
collapse into BD mass objects with circumstellar disks, like more massive proto- 
stars. We rely on numerical experiments with very large resolution, achieved with 
adaptive mesh refinement (AMR). The turbulence simulation presented here is 
the first AMR turbulence experiment ever attempted and achieves an effective 
resolution of 1024 computational zones. The star formation simulation achieves 
an effective resolution of (10 6 ) 3 computational zones, from a cloud size of 5 pc to 



protostellar disks resolved down to 1 AU. 



1. Introduction 

Brown Dwarfs (BDs) are approximately 
as abundant as solar mass stars (e.g. 
iBeiar et all200lHChabrieif2 002). The typ- 
ical Jeans' mass in star-forming clouds is 
of order a solar mass, approximately two 
orders of magnitude more massive than a 
BD. This is usually considered the funda- 
mental reason why solar mass stars are so 
common, but the existence of BDs is a chal- 
lenge to this simple explanation. The origin 
of BDs is apparently an important test for 
the theory of star formation. 



From a different perspective, we could 
instead assume that the origin of BDs is 
quite different from that of more massive 
stars, in which case their existence would 
not necessarily provide any constraints to 
the theory of star formation. This approach 
is exemplified in th e recent suggestion by 
iRemurth Clarke! <j2001l) that BDs are the 
result of the gravitational fragmentation of 
protostellar disks, but that work does not 
draw a definite line separating BDs from 
more massive stars. At what mass does star 
formation switch from the standard mecha- 
nism to the disk fragmentation mode? Most 
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likely not exactly at the hydrogen burning 
limit of 0.075 m©, as this limit is defined by 
nuclear physics and can hardly affect the 
fragmentation process. 

There is also no indication of a discon- 
tinuity in the stellar initial mass function 
(IMF) at approximately 0.075 m©, and no 
evidence of special properties of young BDs 
relative to more massive stars, suggesting 
that perhaps a unique process can explain 
the formation of bot h hydrogen burning 
stars an d BDs (e.g. iJavawardhana et al.l 
2003aE). 

The possibility that BDs and more 
massive stars have a common origin 
due to the process of turbulent frag- 
me ntation has been r e cently proposed 
by iPadoan fc Nordlundl l|2004ft . In that 
work, we show that BDs are predicted 
by our model for the origin of the stel- 
lar IMF from turbulent fragmentation 
llPadoan fc Nordmndl2 002). roughly as fre- 
quently as inferred from observations. Here 
we briefly summarize our theoretical pre- 
dictions and present new numerical exper- 
iments. 

2. Stellar and BD Masses from 
Turbulent Fragmentation 

The gas density and velocity fields in star- 
forming clouds are highly non-linear due to 
the presence of supersonic turbulence. The 
kinetic energy of the turbulence is typically 
100 times larger than the gas thermal en- 
ergy on the scale of a few pc (the typical 
rms Mach number is of order 10) and the 
gas is roughly isothermal, so that very large 
compressions due to a complex network 
of interacting shocks cannot be avoided. 
Under such conditions the concept of grav- 
itational instability, based on a comparison 
between gravitational and thermal energies 
alone in a system with mild perturbations, 
does not apply. Turbulent density peaks 
of any size can be formed in the turbu- 
lent flow, independent of the Jeans' mass. 
Peaks that are massive and dense enough 
(more massive their own Jeans' mass) col- 
lapse into protostars, while smaller subcrit- 



ical ones re-expand into the turbulent flow. 
This is a process that we call turbulent frag- 
mentation 1 , to stress the point that stars 
and BDs, formed in supersonically turbu- 
lent clouds, are not primarily the result of 
gravitational fragmentation. 

Because density peaks formed by the 
turbulent flow need to be larger than their 
critical mass to collapse, a necessary con- 
dition for the formation of BDs by super- 
sonic turbulence is the existence of a fi- 
nite mass fraction, in the turbulent flow, 
with density at least as high as the crit- 
ical one for the collapse of a BD mass 
density peak. As the PDF of gas den- 
sity in supersonic turbulence is fully de- 
termined by the mean density and the 
rms Mach number of the turbulent flow 
iNordlund fe Padoanl Il999t lOstriker et all 
11999ft . we can compute this necessary con- 
dition for BD formation as an integral over 
the PDF of gas density. The r esult, dis- 
cussed in IPadoan fc Nordlundl l|2004f) . is 
that in typical molecular clouds approxi- 
mately 1% of the total mass is available for 
the formation of BDs. This fraction is ap- 
proximately 10% for conditions found in re- 
gions of cluster formation. Based on this re- 
sult, the observed frequency of BDs is con- 
sistent with the distribution of turbulent 
pressure, suggesting turbulence may play 
an important role in the the origin of BDs. 

In IPadoan fc Nordlundl {2002) the stel- 
lar IMF is interpreted as the result of the 
turbulent fragmentation of star-forming 
clouds. It is shown that: i) The power 
law slope, s, of the stellar IMF is directly 
related to the turbulent power spectrum 
slope, (3, s — 3/(4-/3); ii) the abun- 
dance of stars of mass smaller than the 
IMF peak depends on the probability den- 
sity function (PDF) of the gas density in 
the turbulent flow. As a result, the IMF is 
a pow er law with slo pe close to Salpeter's 
value l)Salpeterlll955h at large masses, and 
the abundance of low mass stars and BDs 
is a function of the turbulent rms veloc- 



This term was probably first introduced by 
iKolesnik fc Ogul'Chanskil (|l99Cfl 
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Fig. 1. Left panel: Analytical mass distributions computed for (n) = 10 4 cm -3 , T = 10 K 
and for three values of the sonic rms Mach number, Ms = 5, 10 and 20 (solid lines). The 
dotted lines show the mass distribution for T — 10 K, Ms = 10 and (n) = 5 x 10 3 cm -3 
(lower plot) and (n) — 2 x 10 4 cm" 3 (upper plo t). Right panel: IMF of the cluster IC 
348 in Perseus obtained bv lLuhman et all l)2003[) (solid line histogram) and theoretical 
IMF computed for (n) = 5 x 10 4 cm" 3 , T = 10 K and M s = 7 (dashed line). 



ity, the mean gas density and the mean 
temperature. The IMF is fully determined 
once these physical parameters of the star- 
forming region are known. As shown in 
Figure n the predicted BD abundance in- 
creases with increasing Mach number and 
density (left panel) and is consistent with 
the abundance inferred in stellar clusters 
(right panel). 

3. AMR Experiments 

We employ a 3-D parallel structured 
adaptive mesh refinement (AMR) code, 
Enzo, 2 developed at the Laboratory 
for Computational Astrophysics by 
Bryan, Norman and collaborators 
(iBrvan fc NormarJ Il999t lO'Shea et all 
12004ft . Enzo is a public domain Eulerian 
grid-based hybrid code which includes 
hydrodynamic and N-body solvers and 
uses the AMR algorithm of B erger & 
Colella (|Berger fc Colellal Il989h to im- 
prove spatial resolution in regions where 
that is required. The central idea behind 
AMR is to solve the conservation laws 
on a grid, adding finer meshes in regions 
that require enhanced resolution. AMR 

2 See http://cosmos.ucsd.edu/enzo/ 



is spatially- and time-adaptive, can be 
used with accurate methods for solving 
the (magneto)hydrodynamic equations. 
Mesh refinement can be automatically 
advanced to an arbitrary level, based on 
any combination of problem-specific local 
refinement criteria (e.g., density threshold, 
strong shocks, high vorticity, large grad i- 
ents, Jeans length l)Truelove et al.l [19971. 
cooling time, etc.). 

With this code, we have run the 
first AMR turbulence experiment ever at- 
tempted, achieving an effective resolution 
of 1024 3 computational zones (Kritsuk et 
al. in preparation). The rms Mach num- 
ber of the turbulent flow is M s w 6, and 
we have adopted an isothermal equation of 
state, periodic boundary conditions, large 
scale random forcing and mesh refinement 
based on the local shock or shear ampli- 
tude. The logarithm of the projected den- 
sity field from a snapshot of this simu- 
lation is shown in Figure El The density 
field of this simulation shows an incredibly 
rich structure, with sharp filaments on all 
scales, broken into even denser peaks. Some 
of these density peaks are gravitationally 
unstable and would collapse if gravity were 
included in the experiment. We refer to 
these density peaks as turbulent seeds, as 
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Fig. 2. Logarithm of the projected density field from an AMR simulation of isothermal, 
supersonic hydrodynamic turbulence with an effective resolution of 1024 3 computational 
zones (Kritsuk, Padoan & Norman, in preparation). 



they are the turbulent progenitors of proto- 
stars. Turbulent seeds of BD mass are com- 
monly found in the turbulent flow. 

In order to verify the outcome of 
the gravitational collapse of the turbulent 
seeds, we run a second experiment with 
the same code, where self-gravity is in- 
cluded. In this experiment we use the AMR 
method to resolve the gravitational col- 
lapse. The computational box represents a 
star-forming cloud of size 5 pc and aver- 



age number density 500 cm~ 3 . In the col- 
lapsing regions we are able to fully resolve 
protostellar disks down to a scale of ap- 
proximately 1 AU. This simulation shows 
that the densest turbulent seeds, found 
primarily within dense filaments, collapse 
into protostars surrounded by protostellar 
disks. Figure shows a typical filament 
containing many collapsing "seeds" (left 
panel). The seeds appear as small dense 
cores, because of the limited dynamical 



Padoan et al.: Brown Dwarfs from Turbulent Fragmentation 



5 




Fig. 3. Left panel: Logarithm of the projected gas density in our AMR simulation of star 
formation in a turbulent cloud (see text). Right panel: Higher resolution view of a small 
region inside the dense filament marked by the small square in the left panel. The small 
square is drawn larger than its actual size, as the magnification factor between the right 
and the left panels is 840. 



rage of the image. In fact, visualized with 
higher levels of refinement, they are seen to 
be collapsing into protostars with circum- 
stellar disks (right panel of FigureOJ. Many 
of these protostars contain a fraction of a 
solar mass, some only a BD mass. 

Further analysis and runs with longer 
integration time will allow to establish the 
protostellar mass distribution. However, 
these simulations already teach us an im- 
portant lesson: Turbulent seeds of BD 
mass are common and may collapse into 
BDs, with no need to form BDs from 
the fragmentation of protostellar disks. 
Simulations that do not resolve well the 
turbulent flow (like perhaps every smooth 
particle hydrodynamic (SPH) simulation 
to date) are bound to artificially sup- 
press this important mechanism of BD for- 
mation. Furthermore, insufficient numer- 
ical resolution may also cause artificial 
fragmentation of protostellar disks (Fisher 
et al., in preparation). Recent attempts 
to support the disk fragmentation ori- 
gin of BDs, and even deny the impor- 



tance of turbulent fragmentation, b ased on 
SPH simulations l)Bate et alJl2002j) are af- 
fected by both problems. We also warn 
against a direct application of results 
from simulations without magnetic fields to 
present day star formation, because mag- 
netic fi elds affect the fragmentati on mech- 
anism l|Padoan & NordlundlH999(l and the 
angular momentum transport. 



4. Conclusions 

High resolution numerical experiments of 
supersonic turbulence and of star formation 
in turbulent clouds support the suggestion 
that BDs are formed primarily from the 
process of turbulent fragmentation, as tur- 
bulent seeds, like hydrogen burning stars. 
As frequently pointed out at this work- 
shop, all the observational evidence points 
to a common origin of stars and BDs. 
Furthermore, observations have shown 
that the mass distribution of prestel- 
lar condensations is indistinguishable 
from the stellar IMF l|Motte et alJ Il998t 
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2002), both in the functional shape and 



in the range of masses, i nclud ing BD 
mass cores llWalsh et alJ l200l . This 
direct evidence of the turbulent seeds 
provides strong support to the turbulent 
fragmentation origin of BDs. 
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